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Abstract 



In this paper, we study Raychaudhuri's equation in the background of it! + j3R 2 gravity 
with a phenomenological matter (p oc a(t)~ n ). We conclude that even though the Strong 
Energy Condition (S.E.C.) for Einstein's gravity, which guarantees singularity, is n > 2 for 
p oc a(t)~ n ) a perturbative analysis of Raychaudhuri's equation in the background of R + (3R 2 
gravity reveals that the big bang singularity may not be guaranteed for n > 4. We derive 
the following Strong Energy Conditions for R + (3R 2 ((3^0): 

1) For k < FRW metric, S.E.C. is (0 < n < 4) i.e., -p n < p n < \p n . 

2) For k = FRW metric, S.E.C. is ( 1 < n < 4) i.e., -\p n < p n < \p n . 

3) For k > FRW metric, S.E.C. is (2 < n < 4) i.e., -\p n < p n < \p n . 



I. INTRODUCTION 

An essential mathematical criterion in the singularity theorem is the Strong Energy 
Condition (S.E.C.). The theorem predicts, modulo plausible assumptions, in the existence 
of black holes and cosmological singularity [1,2] 

Singularity theorem is limited by its classical aspect. The question of validity of this 
theorem or how the theorem should be altered when quantum effects are incorporated is a 
fascinating question. 

A straightforward method of incorporating quantum effects is to start from the classical 
Einstein equations and quantize both gravity and matter via corresponding principle, i.e., 
{x,p} P oisson [x,p\qm- Specifically, utilizing Schrodinger representation for canonical mo- 
menta results in the Wheeler-DeWitt equation [3,4]. As for the question of quantum effects 
on the big bang singularity, study of the Wheeler-DeWitt equation has taught us the im- 
portance of another parameter, i.e., boundary condition (initial condition) for the Wheeler 
DeWitt wavefunction[5-8]. The reasoning is quite simple. Quantum mechanical evolution of 
a system is inherently a description of a diffusive system. As a diffusive system, one can only 
hope to evolve the system forward in time from some specified initial (big bang) condition 
on a Wheeler DeWitt wavefunction. There are also models for which the resulting Wheeler 
DeWitt equations become eigenvalue problems and the arbitrariness of initial conditions is 
removed [9-11]. 

Recently, additional method of incorporating quantum effects has come into focus. Ad- 
vances in quantum gravity have revealed that, even if one starts with the Einstein's gravity 
coupled to matter, quantum loop effects of gravity+matter and renormalization procedures 
result in a quadratic gravity [12, 13]. Therefore, even a classical analysis of quadratic gravity 
is inherently a semiclassical analysis. 
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The two methods of incorporating quantum effects, which differ by the choice of canonical 
or covariant quantization, are not necessarily equivalent. This is because a formal equivalent 
of canonical and covariant quantization is only true for renormalizable theories. 

Several authors have studied classical solutions of R + f3R 2 gravity without matter and 
have concluded that Big Bang singularity may be avoided [14-16]. There are several prob- 
lems with higher derivative theories [17,18], e.g., need for additional initial conditions in a 
formulation of a cauchy problem, existence of run away solutions, and whether solutions 
obtained from R + f3R 2 gravity reduce to solutions of Einstein's gravity as f3 — > 0. In order 
to resolve these issues, the author has recently proposed an alternate approach to R + f3R 2 
gravity in which (3R 2 is treated as a backreaction on Einstein's gravity [19]. In essence, 
in this approach the gravitational degree of freedom is not altered from those of Einstein's 
gravity. 

Using such an interpretation of R + f3R 2 gravity, the author also investigated a classical 
and Wheeler DeWitt evolution of it! + f3R 2 gravity for a particular sign of f3, corresponding 
to non- tachyon case. Matter was described by a phenomenological p oc a{t)~ n . It was 
concluded that both the Friedmann potential U(a) (a 2 + 2U(a) = 0) and the Wheeler 
DeWitt potential W(a) ( — + 2W(a) ip{a) = 0) develop repulsive barriers near a for 
n > 4 (i.e., p > |p). The interpretations was clear. Repulsive barrier in U(a) implies that a 
contracting FRW universe (k > 0, k — 0, k < 0) will bounce to an expansion phase without 
a total gravitational collapse. Repulsive barrier in W(a) means that a ~ is a classically 
forbidden region. Therefore, probability of finding a universe with the big bang singularity 
(a = ) is exponentially suppressed. 

Superficially, the prediction of no cosmological singularity for n > 4 (i.e., p > |p) seems to 
be in violation of the Singularity Theorem for Einstein's gravity, which predicts singularity 



for a matter satisfying the S.E.C., i.e., p > — |p and p > —p or equivalently n > 2 for 
p oc a(t)~ n [20]. 

In this paper, we study Raychaudhuri's equation in the background of R + f3R 2 gravity 
coupled to matter (p oc a(t)~ n ). We conclude that the appropriate S.E.C. for R + (3R 2 
gravity is different from that of Einstein's gravity. We also derive explicit expressions for 
S.E.C. for R + (3R 2 gravity. The S.E.C. for R + (3R 2 gravity is shown to be in agreement 
with the author's previous work, which demonstrated that both classical and Wheeler DeWitt 
solutions of R + f3R 2 gravity were free of cosmological singularity for n > 4. 

Sign conventions used in this paper are as follows, g = (—,+,+,+), R ab — \g a bR = 
(+) 8nGT ab , (+) R(U, V) = V V V V - V V V V - V [uy] . 

II. REVIEW OF STRONG ENERGY CONDITION FOR 
EINSTEIN'S GRAVITY 

Even though the goal of this paper is to study Strong Energy Condition (S.E.C.) for a 
quadratic gravity, much of the technique will be borrowed from the analysis of S.E.C. for 
Einstein's gravity. Therefore, we shall briefly review the derivation of S.E.C. for Einstein's 
gravity. For a pedagogical review, please see [20]. 

Let £ a be a tangent vector for a congruence of timelike geodesies. For a hypersurface 
orthogonal congruence, Raychaudhuri's equation is 

^ = -le 2 -a ab a« b + R ab ee. (2.1) 
6 and a ab , are the expansion and sheer of two nearby tangent vectors, respectively. In 
Einstein's gravity, S.E.C. and the ensuing singularity theorem, follow from requiring that 



- R ab Ce = 8ttG 



Tab — -Tg a b 



> (2.2) 



for all timelike £ a . In such a case, ^ < and a pair of nearby timelike geodesic vectors 



converges and will eventually intersect. 

We will be interested in the cosmological singularity. In a FRW metric, matter is de- 
scribed by 



rpa 

2 b — 



( -p 



V 



V 



V 



\ 



V I 



= pt a t b + px a x b + py a y b + px a z b . 



(2.3) 



{t a , x a , y a , z a } are eigenvectors of T b a . They are normalized as — t a t a = x a x a = y a y a = z a z a = 
1. 

Because of isotropy, one can always rotate the coordinate such that the most general forward 
timelike vector is 



£ a = At a +(A 2 -l) 1/2 x a , (AG (l,oo), r& = -l)- 



(2.4) 



Then, 



Rat,?? = SttG 



Tab — -^Tg ab 



a b = 8ttG 



p A 



+ P[A 



(2.5) 



S.E.C. is a requirement on an equation of state (p, p) for which —R ab ^ a ^ b > is true for all 
timelike vectors £°, i.e., A e (1, oo). 

Expression (2.5) is monotonic function of A 2 . Therefore, if we can find an equation of 
state (p, p) for which —R ab £ a £ b > is satisfied by two extreme timelike vectors (A = 1 
and A — > oo), then this (p, p) is guaranteed to satisfy —R a b£ a £ b > for all timelike vectors 
A e (1, oo). Requiring -R ab £ a £ b > for A = 1 gives p + 3p > 0. Requiring -R a b£ a £ b > 
for A — > oo gives p + p > 0. These are the familiar S.E.C. [20]. 
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III. DERIVATION OF STRONG ENERGY CONDITION FOR 
R + R 2 GRAVITY 

Raychaudhuri's equation describes how a congruence of timelike geodesies deviate from 
one another. Therefore, Raychaudhuri's equation (2.1) is valid even in a background of 
quadratic gravity. As in Einstein's gravity, one can be assured of a convergence of congruence 
of timelike geodesies by requiring —R ab ^ a ^ b > 0. In order to proceed, we need the relevant 
"Einstein" equations for quadratic gravity. A pedagogical derivation is given in Appendix. 

X -Rg ah - R ab + 16ttG(3 {^R 2 g ab - 2RR ab + 2R%g ab - 2R, a ,^j = 8nGT ab . (3.1) 
The trace of this equation is 

6 • 16nG(3R% + R = MGT. (3.2) 

Some comments are in order. First, by dimensional consideration, f3 is dimensionless. 
Second, as noted by various authors [21,22], (3.2) resembles a scalar field equation with 
m 2 = —(6 • IGirGpy 1 . Therefore, (3 < is needed to eliminated tachyons. Third, an 
order of magnitude estimate reveals that contributions from the quadratic terms are smaller 
than those of Einstein's terms by (3GR (3G 2 p (3p/p P ianck ■ Therefore, if we make a 
reasonable assumption that (5 ~ 0(1), than we will be justified in treating the quadratic 
terms perturbatively until the very early universe. For a more elaborate discussion on the 
justification for treating the f3R 2 as a perturbation, please see [19]. The subsequent analysis 
should be viewed as a perturbative probe into a possible nonlinear phenomenon. On the 
other hand, it is worth noting that Mijic[23], Starobinsky[24], and Berkin[25] have studied 
the large (3 range and concluded that even for a vacuum, gravity alone can generate an 
inflationary phase in R + f3R 2 gravity. 
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For use in Raychaudhuri's equation, we need an expression for Ricci tensor R ab in terms 
of T ab . 

Combining (3.1) and (3.2) we get 

- R ab + \QixG(3 {^R 2 g ab - 2RR ab - R%g ab - R, a , b ^j = 8uG{T ab - ^Tg ab ). (3.3) 

We are interested in the first order contribution from PR 2 to Raychaudhuri's equation. 
Therefore, for terms already multiplied by (3 in (3.1), we may substituting R ab = —8nG (T ab — \Tg ab j-\- 
0(13) and R = 8nGT + 0{(3) to get 

{^R 2 g ab - 2RR ab - 2R%g ab - 2R, a . b ^j = -G {^T 2 g ab G - 2TT ab G + T%g ab + 2T ;a;fe ) +0{(3). 

(3.4) 

We have introduced a notation G = 8nG. And finally, for use in Raychaudhuri's equation, 
we get 

_ d tatb = G ( Tab ~ \ T 9ab) + 

Kabi i - ^ 2p& (i T 2 gab G _ 2TT ab G + T}™g ab + 2T ]a;b ) 

The rest of the procedure more of less mimicks Einstein's case. Expression for a general 
timelike vector £ a is (2.4). T ab appropriate for a FRW metric is (2.2). 

We will assume that during any epoch in the evolution of universe, the universe is domi- 
nated by a matter with a characteristic dependence on the scale factor (i.e., p — A local 
conservation of T ab gives p = 12 -^-^k ■ 

A purist may argue that a proper way of including matter in the early universe is by 
incorporating a quantum field (e.g., scalar field). We feel that in the study of Raychaudhuri's 
equation, describing matter by a phenomenological p = ^ is adequate. This is because 
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Ce + 0([3 2 ). (3.5) 



Raychaudhuri's equation studies geometrical optics limit of a field in a background metric 
and matter. If we had used a scalar field, than we would have had a difficult task of splitting a 
field (f>(x, t) = (j>o(x, t) + 5(f)(x, t), where (f>o(x, i) and S(j)(x, t) are a background low frequency 
and a fluctuating high frequency components, respectively. The background component 
4>o(x,t) would again result in an effective p = (T® [(/>o(x, t)]) space oc a~ n for some (n). The 
value of (n) would depend on properties of the field (e.g., mass, self coupling V(<f>)). And 
Raychaudhuri's equation would correspond to a geometrical optics limit of 5(f)(x,t). 
Continuing, it is a laborious exercise to show that (3.5) becomes 



n-2. 



G Pn + [3G 2 \3n(n - l)(n - 4)Gp 2 n - 6n 2 (n - 4)ka~ 2 p n 



+ 



(A 2 - I) [lG Pn + (3G 2 \2n 2 (n - 4)Gp 2 n - 4n(n + 2)(n - 4)k a - 2 p n 
A reader who may want to derive (3.6) from (3.5) will find the following helpful. 



(3.6) 



T = T ab g ab = 3p-p=(n-4)p n 



(3.7) 



T ab Ce = A 2 (p + p)-p (3.8) 
T, a ,re = A 2 (d 2 T - -d t T) + -d t T (3.9) 

T% = -a- 3 d t (a 3 d t T) = -(n - A) a - 3 d t (a 3 d t p n ) (3.10) 

Again, Strong Energy Condition (S.E.C.) is requirement on equation of state for which 
—R a b£, a £, b > for all forward timelike vector £ a . There is a subtlety to keep in mind. From 
(2.1), even though —R a b(, a t; b > implies that a pair of timelike geodesies converge, the 
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opposite —R a bi a i h < does not imply that a pair of timelike geodesies diverge. This is 
because of existence of other negative negative terms in (2.1) 

Expression (3.6) is a monotonic function of A 2 . Similar to Einstein's case, if we can find 
an equation of state (p, p) for which —R ab ^ a ^ b > is satisfied by two extreme timelike vectors 
(A = 1 and A — > oo ), then this (p, p) is guaranteed to satisfy —R a bi a i h > for all timelike 
vectors A e (1, oo). 

Requiring —R a bi a i b > for A = 1 gives 



n - 2 ~ 



+ 3n(n- l)(n-A)Gp n - 6n 2 (n - 4)A;a _ > n > 0. (3.11) 



2 

Requiring —R ab £ a £ b > 0. for A — > oo gives 

-Gp n + /3G 2 2n 2 (n-A)Gp 2 n -An{n + 2)(n-A)ka- 2 p n > 0. (3.12) 
3 

The problem of cosmological S.E.C. has been reduced to solving for (n) that will satisfy 
(3.11) and (3.12) as a -> 0. 

As a partial check in algebra, consider the case of f3 — 0. (3.11) gives n > 2 , which 
is equivalent to p n > — |p n . And (3.12) gives n > 0, which is equivalent to p n > — p n . As 
expected, these are the S.E.C. for Einstein's gravity. 

Now for p ^ case (f3 < for a tachyon free system), the argument is significantly more 
subtle. We note that as one goes further back in time, depending on the value of (n), one of 
the (p n , p^, a~ 2 p n ) in (3.11) and (3.12) will grow fastest and hence dominate the expressions. 

We already know the result if p n dominates (i.e., (3 = 0). We need n > 2 to satisfy (3.11) 
and (3.12) (Figure la). 

Now if p 2 n dominates, than (3.12) reduces to f3n(n — l)(n — 4) > 0, i.e., 1 < n < 4 since 
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P < . And (3.12) reduces to f3n 2 (n — 4) > 0, i.e., < n < 4. The two are simultaneously 
satisfied by 1 < n < 4 (Figure lb). 

And finally, if the curvature term {ka~ 2 p n ) dominates, than (3.11) reduces to —kf3n 2 {n — 
4) > and (3.12) reduces to —k(3n(n + 2)(n — 4) > 0. The two reduced equations are 
satisfied simultaneously for (k > 0) FRW metric by n > 4 (Figure lc), and for (k < 0) FRW 
metric by < n < 4 (Figure Id). We have tactfully assumed that n > 0, i.e., energy density 
should not get less dense when squeezed. 

We also note that only when < n < 2 can the curvature term ka~ 2 p n grow to dominate 
over the p 2 n term, as one goes further back in time (a — > 0) 

We have obtained various collection of results. A slightly more useful conclusion would 
be, for a given (p n ,k,f3) whether there will be a big bang singularity. We will proceed to 
address this question in two steps. The strategy will be as follows. First, for a given (p„, k, (3) 
we will have to determine which of the (p n , p 2 w a~ 2 p n ) will dominate (3.11, 3.12)) as a — > 
. Then from Figures lb-Id, we will be able to read off whether such a given set of (p„, k, (3) 
has a big bang singularity. 

Let us first consider k = FRW metric, with p oc a(t)' n . As one goes back to earlier 
time, p 2 n should dominate over p n . Therefore, from (Figure 2a) the appropriate S.E.C. is 
(1 < n < 4). 

A careful combination of previous results reveals that the S.E.C. for k < FRW metric 
is < n < 4 (Figure 2b). The reasoning is as follows. As a — > 0, either ka~ 2 p n or p 2 n could 
end up dominating (3.11, 3.12). For < n < 2, ka~ 2 p n will dominate. And from Figure Id, 
all of < n < 2 results in a convergence of a pair of timelike vectors. For n > 2, p 2 n will end 
up dominating. From Figure lb, of this region, only 2 < n < 4 results in a convergence of a 
pair of timelike vectors Q.E.D. 
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Similarly, analysis reveals that the S.E.C. for k > FRW universe is 2 < n < 4 (Figure 
2b). Again for < n < 2 , ka~ 2 p n will dominate. But from Figure Id, none of this region 
results in a convergence of a pair of timelike vectors. For n > 2, p 2 n will again end up 
dominating. And from Figure lb, of this region, only 2 < n < 4 results in a convergence of 
a pair of timelike vectors Q.E.D. 

These results, of cosmo logical strong energy condition for R + R 2 gravity, are summarized 
in Figures 2a-2c. 

IV. DISCUSSIONS AND CONCLUSION 

We were primarily interested in how Big Bang singularity would be affected by quantum 
effects. Granted that even a classical analysis of R + f3R 2 gravity can be interpreted as 
a semiclassical analysis, there are several limitations to our analysis and we would like to 
briefly raise these points. 

1) Near the Planck epoch, large quantum fluctuations will undoubtedly result in an 
inhomogeneous universe, yet we have assumed homogeneous (and isotropic) metric in (2.3- 
2.4) to make the problem analytically tractable. A future work will have to address how the 
present conclusion is affected by anisotropy and inhomogeneity. 

2) Raychaudhuri's equation is a classical geometric optics limit, which may be invalid 
near the Planck epoch. 

3) In Raychaudhuri's equation, we have taken the limit as a — > 0. A more physical limit 
might be to cut off the limits at p n — > p p i an ck, ka~ 2 — > {planck length)" 2 . Various set 
of conclusions obtained by this limiting process turns out to be sensitive to fine tuning of 
parameters and will not be discussed here. An interested reader is invited to explore various 
possibilities. 
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4) In retrospect, in the deriving of S.E.C. for R+f3R 2 gravity, we have pushed the analysis 
until kf3a~ 2 p n or fip 2 n dominated over p n . Yet, if the terms linear in (5 ever grow larger than 
the terms from Einstein's gravity, than this will be outside the validity of perturbation. This 
is a serious objection. But as it is a common practice in perturbation theory, we are inclined 
to interpret these results obtained by a perturbative analysis as a possible glimpse into a 
nonperturbative effect. With this in mind, we must settle for a following weaker conclusion. 

"Cosmological S.E.C. for Einstein gravity is n > 2 . For R + f3R 2 gravity, first order 
perturbation analysis indicates a deformation of S.E.C. region such that: 

1) for k < 0, n > 4 may have to be excluded and < n < 2 may have to be included. 

2) for k = 0, n > 4 may have to be excluded and 1 < n < 2 may have to be included. 

3) for k > 0, n > 4 may have to be excluded." 

In closing, these results compliment the authors recent work, which demonstrated that 
both classical and Wheeler DeWitt solutions of R+/3R 2 gravity were free of big bang singular 
free for n > 4. 

V. APPENDIX 

The action for quadratic gravity is 

1 = ~T^G I d4x ^~9 R ~ J d ^ x [P lR2 + fo R ^ Rab + foRabcdR abcd ] + Imatter + surface term. 

(5.1) 

We have formally included a surface term to cancel any boundary term that would result 
in applying the variational principle. We will be interested in applying the formalism to 
homogeneous and isotropic metric, i.e., Wely tensor vanishes C a b c d = [26]. By definition of 
Weyl tensor, C a bcdC abcd = R a bcdR abcd — 2R ab R ab + \R 2 . This gives one relationship among 
the possible quadratic terms. 

11 



The second relationship is from the 4 dimensional generalization of Gauss Bonnet formula 
[12], 



R 2 - AR ab R ab + R abcd R abcd = exact derivative. (5.2) 

The two relationships, combined with the fact that Euler Lagrange equations are unchanged 
by addition of an exact differential, allow any two of Pi, @2, flz to be set equal to zero in the 
action (5.1). We choose to set (3 3 = (3 2 — 0. 
The resulting Euler Lagrange equations are 

X -Rg ah - R ab + WttGP {^R 2 g ab - 2RR ab + 2Rf a g ab - 2R, a . b ^j = 8nGT ab . (5.3) 

The trace of this equation is 

6 • I&kGPR* + R = 8ttGT. (5.4) 
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FIGURE CAPTIONS 



Figure a: Strong Energy Condition for Einstein's gravity (i.e., (3 — 0). 
Figure b: —R a b^ a ^ b > determined under the assumption of p\ domination. 
Figure c: —R a bt; a t; b > determined under the assumption of ka~ 2 p n domination 
(k > 0). 

Figure d: — -R a &£ a £ b > determined under the assumption of ka~ 2 p n domination 
(k < 0). 

Figure a: Strong Energy Condition for k = FRW metric, (1 < n < 4 i.e., — |p n < 
Pn < |Pn )• For 1 < n < 4, a universe ends in domination. 

Figure b: Strong Energy Condition for k < FRW metric, (0 < n < 4 i.e., — p n < 
Pn < |p«)- For < n < 2, a universe ends in ka~ 2 p n domination. For 2 < n < 4, a 
universe ends in p^ domination. 

Figure c: Strong Energy Condition for k > FRW metric, (2 < n < 4 i.e., — |p„ < 
Pn < |Pn)- For 2 < n < 4, a universe ends in p\ domination. 
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